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Abstract. We consider sets of quasimodes for the Dirichlet Laplacian on a domain with boundary 
where the geodesic flow exhibits mixed dynamical behavior. We assume that the billiard flow has 
an invariant ergodic component, U, and study defect measures, \i, of positive density sequences of 
almost orthogonal quasimodes. We demonstrate that these measures have j.t\ u = c^ L \ u where fa 
is the Liouville measure. In order to do this, we adapt quantum ergodicity results to the case of 
quasimodes. Finally, using Bunimovich's mushroom billiards |2|, we give an example where our 
results apply. 



1. Introduction 

The distribution of eigenfunctions over phase space in the semiclassical limit is an important 
object of study in the theory of quantum chaos. The fundamental result is a quantum ergodicity 
theorem of Shnirelman [14], Zelditch |fT5l , and Colin de Verdiere [4j which states that, for clas- 
sically ergodic systems, high energy eigenfunctions distribute uniformly in phase space. Since 
we will study domains with boundary, we note that Zelditch-Zworski lH6ll generalized quantum 
ergodicity to that case. 

Some progress has been made toward understanding semiclassical limits of eigenfunctions 
for systems with dynamical behavior that is not completely ergodic. Marklof and O'Keefe, [8], 
examine separated phase spaces for certain maps. More recently, Marklof and Rudnick [9| have 
made further strides towards an understanding of quantum behavior for mixed phase space. In 
particular, they prove that, for rational polygons, the eigenfunctions of the Dirichlet Laplacian 
equidistribute in configuration space. We use [9| as inspiration for further work on semiclassical 
limits for systems with mixed dynamical behavior. 

We examine systems whose phase spaces have invariant subsets, U, such that the flow restricted 
to U is ergodic and the Liouville measure of U is positive. In particular, let (M,g) be a smooth 
Riemannian manifold of dimension d with a piecewise smooth boundary, dM. Then M c M, 
where M is a manifold without boundary to which g extends smoothly, and dM = Uy =1 Nj where 

Nj are smooth embedded hypersurfaces in M. Define d°M c dM to be the open set of all points 
where the boundary is smooth. Then the complement dM \ d°M has measure zero. 

We consider P(h) — —h 2 A + V with Dirichlet boundary conditions and let p(x, E) be its principal 
symbol. Since p(x, £) is smooth up to the boundary, we can extend it smoothly to T*M. We assume 
that 



(1.1) 



for E € [a,b], dpL-i^ + 
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Figure 1. We show two billiards trajectories in a Bunimovich mushroom [2 J with 
semicircular hat of radius 1 and centered base of width 1/2 and height 1. On the 
left, we have a trajectory in the ergodic portion of phase space. On the right, we 
have one in the integrable portion. 



and 

(1.2) x € d°M, V(x) = E^>dV t N* x d°M. 

so that p -1 (£) and T*„ M M intersect transversally. 
We then write 

p-\E) n T* 0M M = Q+ U Q- U Q°, 

where (x, £) lies in if the vector H p x e TM points outside of M, in Q~ if it points inside M, 
and Q° if it is tangent to dM. The set contains the glancing covectors and, under {L.2) , has 
measure zero in p -1 (£) n T*„ M M. 

We define the broken Hamiltonian flow as follows, (see, for example, [5. Appendix A]). We 
denote this flow by (pt := exp(Hpt). Assume without loss of generality that t > 0. We consider 

exp(tH p )(x, £) defined on T*M, and let fo be the first nonnegative time when exp(tH p )(x, £) hits 
the boundary. If this happens at a non-smooth point of the boundary, or if exp(toH p )(x, e Q°, 
then the flow cannot be extended past i = fo- Otherwise, exp(foH p )(x, E,) e and there exists 
unique (xq, ^o) £ ^£ such that the natural projections of exp(foH p )(x, E,) and (xq, Eo) onto T*dM are 
the same. We then define cpt inductively, by putting (p t (x,E,) = exp(tH p )(x,£,) for < f < fo and 
<pt(x, £) = <p t - to (xo, E,q) for f > fo. Define for T > 0, the set 

S T c rMnp-^fl,!)]) 

to be the closed set of all (x, E) such that <pt(x, £,) intersects a glancing point for some f e [-T, T]. 
Then, as shown in |Tl6l Lemma 1], St has measure zero in p _1 (E). 
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Figure 2. We show two high energy eigenfunctions for the laplacian in a Buni- 
movich mushroom Q. The left hand eigenfunction spreads uniformly through 
the ergodic protion of phase space, and the right hand eigenfunction concentrates 
in the integrable portion. (Images courtesy of A. Barnett and T. Betcke HI). 



The ease with which quasimodes can be constructed in examples motivates us to generalize 
existing quantum ergodicity results to that setting. Although this is a natural generalization 
of quantum ergodicity results, we are not aware of any reference in the literature. We present 
this generalization in section [2] In order to do this, we make the following definition (with the 
obvious analog in the classical setting under the rescaling h = A -1 ). 

Definition 1. A positive density set of quasimodes for P on [a,b], < a <b is a collection {(\pj,Ej),j = 
1, ...} satisfying, 



Remark: Notice that the set of all eigenfunctions is a positive density set of quasimodes. Hence, 
all of our results will apply in this case. 



(i) my = i 



(2) \\(P-E j )4> j \y = o(h 2 < i+1 ) 

(3) \(iPj,ip k )\ = o(h 2d ) j*k 

(4) \{Eje[a r b])\ >ch~ d r c>Q 




(1.3) 



<pt is ergodic on p 1 (E), E e [a, b] 
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we have the following analog of quantum ergodicity for a positive density set of quasimodes 
for P. For a semiclassical pseudodiffential operator of order m, we write B £ W™, and denote its 
semiclassical symbol, 0h(B) (see, for example, 1(171 Section 14.2] for definitions). 

Theorem 1. Suppose that (M, g) is a compact manifold with a piecewise smooth boundary and \1.3) holds. 
Let {(ipj,Ej), j = 1,...} be a positive density set of quasimodes of the Dirichlet realization ofP on \a,b\. 
Then there exist a family of subsets 

A(h) c {Ej\a < Ej < b] 

of full density such that for all B £ ^(M ) with symbol, Oh(B), compactly supported away from dM, 

(Bip jr iPj)^f o h (B)dp Ej 

Jp-l(Ey) 

for Ej £ AQi). 

In section [3J we specialize to the case P - —h 2 A and pass from semiclassical to the classical 
setting. For a classical pseudodifferential operator of order m, we write B e W m and denote its 
classical symbol, o(B) (see, for example, [|6j Chapter 18] for definitions). 

Theorem 2. Suppose that (M,g) is a compact manifold with a piecewise smooth boundary and ( |1.3| 
holds. Let {(ipj,Ej),j = 1, ...} be a positive density set of quasimodes of the Dirichlet realization of -A on 
[0, 1]. Then, there is a full density subsequence ipj k such that for all classical pseudodifferential operators 
B £ V F°(M°) with symbol compactly supported away from dM, 

(Bip jk ,ip jk )^f o(B)dp L . 

JS'M 

In section |3J we continue to restrict our attention to the classical setting. But, instead of the 
ergodic assumption ( |1.3| l, we make the mixed dynamical hypothesis 

(1.4) 311 c S*M, in{U) > 0, ii L {dU \ U) = 

(1.5) (pt is ergodic on U, (pt(U) = U. 

In this case, we have the following result 

Theorem 3. Suppose that (M, g) is a compact manifold with a piecewise smooth boundary and that \1A) , 
\1.5) hold. Let {(ipj, Ej), j = 1, ...} be a positive density set of quasimodes of the Dirichlet realization of -A 
on [0,1]. Then there is a full density subsequence ofipj such that if any further subsequence ipj k has a 
defect measure, \i, p\u = cp.i\u- 

Here, a defect measure for a sequence \pj is the limit of the semiclassical measures pj associated 
to ipj. (see, for example, [17, Chapter 5] for details) In the boundaryless case, Riviere [12] recently 
proved a theorem on accumulation points of semiclassical measures which gives similar results 
for separated phase spaces. He also provides examples of separated phase spaces for the geodesic 
flow on boundaryless manifolds. 

Remark: Hypotheses ( |1.4| include the case where U is closed with positive measure boundary. 
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Finally, in section 5j we apply our results to the Dirichlet Laplacian for mushroom billiards. 
Figure [l] shows two billiards trajectories, one in the ergodic region and one in the integrable re- 
gion. Percival Ml3l made conjectures on eigenfunctions that were numerically verified by Barnett 
and Betcke in [1 J in the case of mushroom billiards. Our results provide rigorous proofs for some 
parts of these conjectures. (For a complete description of the billiard map on mushrooms see [2J.) 
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2. Quantum Ergodicity for Quasimodes 

In this section, we demonstrate how to adapt quantum ergodicity results for eigenfunctions to 
positive density sets of quasimodes. We make the following definition. 

Definition 2. Let E;,j - 1, ... be all the eigenvalues ofP with multiplicity. A complete set of quasimodes 
for P is a collection {ipj,j - 1, ...} satisfying 

\\iPj\y = 1, ||(P - Ej)^j\\ L 2 = 0(^ +1 ), \(lpj, xf, k )\ = o{h 2d ), j * k. 

We will need the following lemmas relating various spectral quantities for operators to their 
corresponding expressions involving complete sets of quasimodes. 

Lemma 1. Let [u^,k = 1,...} be the eigenfunctions ofP corresponding to E^ and let = 1,...} be a 
complete set of quasimodes for P. Then, 

\E r E k \<Ch^ 

Proof. First, observe that, near an energy level E, there exists {c\{h), C2Q1)) with \c\(h) —C2(h)\ > ch d+1 
such that ((E - c 2 (h),E- Ci{h)) U (E + Ci{h),E + c 2 (h))) n {E jr j = 1, ...} = - if this were false, then 
the spectrum of P would violate the Weyl law. 

Now, let A = {ipj\Ej e (E - Ci(h), E + C\{h))) and TI be the spectral projection onto [E - Ci(h), E + 
ct(h)]. By the Weyl law, |A| < Ch~ d . Therefore, by Proposition 32.4], 

dim (IT span (A)) = |A|, 

since 

o{h 2d ) + o{h 2d+l )0{h- d - 1 ) = o(h d ). 

But, rank n = |A|. Hence, span A = range n and the result follows from the almost orthogonality 
of ipj. □ 
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Lemma 2. Let Abe a Hilbert-Schmidt operator and {ipj,j - 1, ...} be a complete set of quasimodes for P 
and a <b. Then, 

(2.1) h d £ IIA^II 2 = h d \\U [aM A\\ 2 HS + o(l) 

a<Ej<b 

and if A is of trace class, 

(2.2) h d Yj (Wp $j) = h d TrU [aM A + o(l). 



a<Ej<b 



Proof. First, let (cj0 = (ufcijjj). Then, by Lemma [TJ 

Jt fc 
where all sums are taken over A: such that — Ej\ < Ch d+1 . Observe that 

k k,k> 
Hence, by the orthogonality of Uj, 

Y c i^' = 8 V + °^)- 

k 

Now, 

h d Y <A$jMj) = h d Y c jk cJ(Au k ,Au' k ) = h d Yffa? + o{h d )){Au h Au v ) 

a<Ej<b j,k,k' k,k' 

= h d \\Tl [aM A\\ 2 HS +h d Y J o(h d )(Au k ,Au k ,) = h d \\U m A\\ 2 HS + 0(1) 

kjc? 

where the last equality follows from the fact that there are at most 0(h~ d ) terms in each sum. 
An analgous argument shows that ( |2.2| holds. □ 



To prove Theorem]!} we need the following lemma similar to [5, Lemma A.2] 
Lemma 3. Let x e Q°(M°) with < * < 1. Then for a' <a<b<b', 

(2nh) d Y f (l-x)tyj\ 2 dVol< f (l-^ + O (l)« S ^0. 

cTTTm^M JTMrarHU'VYi 



Ej£[a,b] ' 



Proof. This follows directly from [5, Lemma A.2] and Lemma |2]once we observe that 

(2nh) d Y f (! - xMjfdVol = (2nh) d Y <(1 - xWj, </>/> = (2nh) d Tr U [aM (l - X ) + o(l). 

Eje[a,b] ^ M E } e[aJ>] 

□ 
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The following lemma allows us to take arbitrary collections of orthogonal quasimodes and 
complete them. Hence, we only have to prove results for complete sets of quasimodes and they 
follow for positive density sets. 

Lemma 4. Let {{ipi, Ey), j - 1,..., J} be a set of almost orthogonal quasimodes with 

\\4>j\y = h \\(P-Ej)ip j \\ L2 =o(h 2d+1 ), \{xPj,^ k )\=o(h 7d ), j*k. 

Then, there exists {(pi, k = 1, ...K} such that the set {(p k ,k = 1, ...,K] U {ipj,j - 1, ...,/} is a complete set of 
quasimodes for P. 

Proof. Let Uy be eigenfunctions of P corresponding to Ey. Select an energy E. Then, by the 
proof of Lemma [lj there are gaps at distance c\(h) of size h d+1 in the spectrum of P near E. Let 
A = {u k : |Ey - E| < ci(h)}. Let A' := {ipj : |Ey - E| = o(h d+1 )}. Then, for ipj e A', we have that 

i/>y = ^ c jk u k . 

u k eA 

Now, define N := |A| and M = |A'| letting bj = (cji,Cj2,...,cp), j - 1,..M. Let {e\,...,e m } be an 
orthonormal basis for span ({by : j - 1, ...,M}). Apply the Gram-Schmidt process to obtain an 
orthonormal basis [e\, e^r, Vm+i, Vn} of R N where 

v k - (v k i,...,v kN ). 

Then, letting 

(p k = Vk i u i' M + l <k <N, 
i 

{ipi, ipMttyM+i, ■■■,<Pn} is an almost orthonormal basis for span A. Repeating this process for 
each cluster, we obtain a complete set of quasimodes. □ 

We also need the following restatement of results in [3. Section 4.3] that is found in O Lemma 
Al]. 

Lemma 5. Fix T > 0. Assume that A e ^ com P(M°) is supported away from the boundary of M and 
WF h (A) c p~ l {[a, b])\S T . Then, for each x £ Q°(M°) and for each t € [-T, T], the operator xe~ itp/h A is a 
Fourier integral operator supposted away from dM and associated to the restriction ofcpt to a neighborhood 
ofWF k (A) n <p~ l (supp x)i V^ us an Ol 2 ^l 2 (^°°) remainder. The following version ofEgorov's Theorem 
holds: 

xe up/h Ae -up/h = Atx + o(^) L2(MHL2(M) , 
where Aj iX e W com P(M°) is supported away from dM and o k (Ar iX ) = xi ao( Pt)- 
Now, we prove Theorem [TJ following [5, Appendix A] 

Proof. Let (\pj, Ey) be a positive density set of quasimodes. By Lemma |4l and the fact that our 
original quasimodes have positive density, we may assume without loss that (ipj, Ey) is a complete 
set of quasimodes (see definition [2]). 
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We first show that 



lim sup h d y 

Ej£[a,b] 



(Btp jt xl>j) - f o h (B)d^ 

Jp-HEj) 



= 0. 



Then, by a standard diagonal argument that can be found, for example, in [[171 Theorem 15.5], 
we extract the set A(h). 



Take a', V such that a' < a < b < V and the assumptions ( |1.2| | and ( |1.1[ ) hold for E 6 \a' ,b'\. 
(If they do not hold when E [a, b], then we need to take a' , V getting close to a and b. e.g. Let 
a' — a — 1/T and b' = b + 1/T. Then estimate the extra contributions by the Weyl law.) Now, fix a 
T > and choose xt e Q°(M°) with < xt < 1 and 



X 



rMnp-^a'-U'-i]) 



(l-rf^a^T- 1 . 



Let i/^eQV -1,&' + 1) have 



i/>(E) f Xrdii E = \ o h (B)d{i E , Ee[a',b']. 

Jp-HE) Jp'HE) 

Then, by Lemma |3j it is enough to show that the conclusion holds for B - ip(P(h))xT, whose 
symbol integrates to on p _1 (E). Therefore, we assume, without loss, that 



X 



G h (B)d[i E = 0, Ee[a',b']. 



p-i(E) 

By elliptic estimates, we may assume that WF;,(B) c p _1 ((a',&')). More specifically, B e y$/ com P. 
Thus, we can write B = B' T + Bj, where WF h (B' T ) Pi S T = and \\o h {B^)\\ L 2 {p -i {a , ih , }) < T _1 . 

Then, by Holder's inequality, Lemma |2]and [5. Lemma 2.2], we have that 

( I W 2 



h d \{B'^j^j)\<C h d 

EjE[a,b] y EjE[a,b] 



c 



+ 0(1) 



h d Y W B T U i\ 

< c\\0 h (B'i)\y {rHa , /hl]) + o(i). 



Bj. 



Hence, the contribution of Bj goes to in the limit limj^oo limsup^Q and we may replace B by 
Now, by Duhamel's formula and the unitarity of e ltp l h , 

"Jo 



So, defining 



(Ah 



■- I f e itP/h Ae -itP/h dt 
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and using Lemma [2] and Holder's inequality, we have 
h d £ KB^V;)! = h d £ \(B' T e- itE i%,e- itE i%)\ = h d £ \((B' T ) T ^p ^)| + o T (h 2d ) 



Eje[a,b] 



Eje[a,b] 



Eje[a,b] 



< 



h d Yj (\\(B' T )Tipj\\+o T (h 2d )) 2 



Nl/2 



E,6[fl,fc] 



xl/2 



Eye[fl,fo] 



+ o T (l) 



From this point forward, the proof is identical to that in |5j Theorem A.2]. By Lemma|5j (B' t )tXt 
is, up to an 0(h°°) L 2^ L 2 remainder, a pseudodifferential operator in vj/ com P compactly supported 
inside M° and with principal symbol 



o h ((B' T hXT) 



XT 



I a h (B' T )o(p t dt. 
Jo 



Now, all that remains to show is 

lim limsup^ Y \\(B') T Uj\\ 2 = 0. 

Eje[a,b] 

Since, by Lemma [3j 

limsup(27i^ Y \\(1-XT)uj\\ 2 L2 <T-\ 

h ~^° Eje[a,b] 

we can replace (B' t )t by (B' t )xt- Thus, by [5, Lemma 2.2], it remains to show that 
( 2 ' 3 ) \\0h((B' T )TXT)\y(p-\[a',b'])) ^ j^m IKo-/ ! (B / T ))rllL2(p-i([a',6'])) = 0. 

To do this, write 

IK^(B / t ))tIIl2( p -i(e)) ^ ll<fffc(B)>rlb(p-i(E)) + W(oh(Bj)h\\m P -HE)y 

The first term on the right goes to when T — > co by the von Neumann ergodic theorem and the 
second term is bounded by ||a^(B^)|| L 2( p -i( E )) and hence also goes to 0. 



□ 



Remark: Notice that ( |2.3| l is the only step in which the ergodicity of the flow is used. This will 
be important when we adapt the result to ergodic invariant subsets of phase space. 

3. From Semiclassical to Classical Quantum Ergodicity 

For completeness and to present the proof in the quasimode case, we will now pass from 
Theorem QJwith P — —h 2 A to Theorem0 

Proof. Let h 2 A 2 = E, and A = h~\ x e C°° X (Q = for |£| < 1 and x = 1 for \E,\ > 2, and Xe = *(£/e). 

Let A be a classical pseudodifferential operator of order on M with symbol compactly 
supported away from dM. Define :— o(A). Then, a$ is homogeneous degree on T*M. Hence 
a (x, D)x(D) = a (x, hD)x(hD/h). Now, define A e e W°(M°) by 

A e := a(x,hD)Xe(hD). 
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Theorem [I] gives, for < a < b, 

h d V (Ae^^-f OhiAeWE, 



0, ft = A" 1 -> 0. 



hAj£{a,b] I J P ( E /) 

But, since Oh{A e ) = ao(*/ £)/£e(£) and flo is homogeneous degree 0, we have 



f o h (A £ )dp E = f a(A)d^ L + 0(e). 

Jp-HEj) JS'M 



Hence, we need to show, 



limlimsupfr d V" ((A - A e )^j, ii>/) = 0. 



fr->0 



hAj€[a,b] 



By Holder's inequality, Lemma |2j and © Lemma 2.2], 

f / 



^ £ - A e )i/>y, \pj)\ < C h d IK^-^/l 

hAj£[a,b] y hAj£[a,b] 

< C\\o h (A - A e )\\ L 2 {p -i ({am + o(l). 

But, linv-»o \\oh(A - A e )\\ L 2 ir i {[aM)) = 0. 
Thus, for any 5 > 0, we have 



C 



\l/2 



hAj£[a,b] 



+ 0(1) 



lims 



supft rf V <A> ; vi/>/>-f 

-0 hAy^l] Js * M 



= 0. 



0. 



All that remains is to show that 

limlimsup^ V (Aipj,ipj)-f a(A)d[Xi 

But, letting Hq - -f a{A)dp. E , and applying the Weyl law, and Holder's inequality, 

\l/2 



h d Yj K%^;>-«ol < C h d Y IK^-W; 
ftA ; e[0,6) ^ fcA ; e[0,6) 



O(6 d )||(A-fl- )|| L2 ^ L2 + O 6 (/ 2 ). 



Once again, to obtain the full density subsequence, we employ a standard diagonal argument. 

□ 

4. Quantum Ergodicity for Subsets 

Now, we prove Theorem [3] using the fact that only step ( |2.3| uses the ergodicity of the flow. 
For completeness, we include the following elementary lemma. 



Lemma 6. Let U satisfy ( |1.4[ ), let \ii be a finite measure on S*M and suppose that for all a e C°°(S*M) 
compactly supported away from dM with J u adp.i = 0, we have J u adfi2 - 0. Then, p2\u - cp.L\ufor 
some c > 0. 
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Proof. Let {i\ = Then, y.\(U) - 1. Define x e C°°(S*M) compactly supported away 

from dM with < x ^ 2 and = 1 and xdy-i = ci > 0. (To see that such functions exist 

simply take non-negative approximations to \\j with support inside M° '.) 

Now, let a e C°°(S*M) compactly supported away from dM with J u adp-i = a. Then, f u a - 

axdf-ii - and hence f u a - ax&\ii - 0. Therefore, for a e C°°(S*M) compactly supported away 
from dM, 

Ju Ju 

But, /Jilu and /^lu are positive distributions of order since U satisfies ( |1.4| . Hence, since 

I fld/j2 = I fld(c2/Ji) = 
Ju Ju 

for all fl £ C°°(S*M) compactly supported away from dM, and M can be exhausted by compact 
subsets, /i 2 |u = c 2 y-i\u- □ 

We now prove Theorem [3] 

Proo/. Lemma |4| shows that without loss of generality we may work with \pj forming a complete 
set of quasimodes. Fix a € C°°(S*M) compactly supported away from dM with 



(4.1) 



I ad\ii 
Ju 



Fix e > 0. Then, let U £ be open with U cU £ and /jl(lI £ \ U) < e (note that we use ( |1.4) here). Let 
^ e £ Q°(ILJ compactly supported away from <?M with Xelu = 1, < x e < 1. Then, let a e = Xe®- 
Based on the arguments used to prove Theorems [T] and |2j we have that 

(4.2) limsupA^ Yj \(AeXpj,ipj)\ < C\\(a e ) T \\ L 2 (s , M) , 



A ->°° Aj<A 



i r T 

«e>T := ~ a e o(p t (x,E,)dt. 
1 Jo 



where A e = a e (x,D), 

i r T 

( 

* JO 

We have that <j9f is ergodic on IT and U is invariant under <p t . Hence, by the von Neumann 
ergodic theorem (a e lu)r — > J u a ed[^L = in L 2 . Again, by the ergodic theorem, (fl e lu c )r — > P«e in 
L 2 with ||Pfl e || L 2 < C||fl||L~e. Hence, 

(4.3) lim |Kfle>Tlb(s-AD ^ C||«||i~e. 

We now show that one can extract a full density subsequence ipj k such that for all a € C°°(S*M) 
satisfying ( |4.1[ ) 

lim lim (A e ip jk ,ip jk ) = 0. 

e— >U /r— >oo 

To do this, fix e > and let 

He) = {Aj < e- 1 : \(A £ ipj,ipj)\ > (C||fl|| L -e)*}. 
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e d \T(e)\ < (Ce||«||i-)$ 



Then, by the Chebyshev inequality, ( |4.2[ ), and ( |4.3[ ), 
and for Ay £ T(e) 

KA^-,^)! <(Ce||fl|| L »)5. 

But, by the Weyl law, 

|r(e)|/|{A ; - < e" 1 }! = e d \T(e)\/ (Vol(ZTM) + o(l)) < C(||a|| L «e) 1/2 . 

Therefore, A(e) = {Aj < e -1 } \ T(e) has full density as e — > 0. Now, we take a countable dense 
set, {%} of C°°(S*M) n {a : ||fl||L~ <1, a satisfies ( |4.1| |} and apply a variant of the standard diagonal 
argument (contained for example, in [17, Theorem 15.5]). 

Now, suppose that ipj k is a further subsequence with defect measure [i. Then, for a e C°°(S*M) 
satisfying ( |4.1[ ), we have that 



k 

Hence, by the dominated convergence theorem, 



c-^oo J 



Jad\i - 0. 
u 



Thus, we may apply Lemma [6] to obtain the result. □ 
Remark: We may also make similar arguments to obtain a result in the semiclassical setting. 

5. An example 

We now present an example to which Theorem [3] applies. Let Q be a symmetric mushroom 
billiard, as in (23, composed of a hat that is a semicircle of radius 1 and a base that has width w 
and height h (figure [l] shows such a mushroom with two billiard flows). Then S*Q has a subset 
U satisfying ([L4) and ([L5]l (see 0). Let 

u nk = sin(n6)J n (a n/k r) 

be the eigenfunctions of the Dirichlet Laplacian on the semidisk. Here } n denotes the Bessel 
function of the first kind of order n and a„ r k the k th positive zero of } n . 

From [10, Appendix A], we have that for < y < |, z 6 (0, 1 - r0'~i) 

< Jn(nz) < 2-" )V3 . 

Also, from [11], we have for < k < c^n, that n < a n ^ < Cn for some C > 0. Thus, } n (^ n ,k r ) = 
0(2-" , /3 )forr < ^. 

Now, suppose < w < Then, let x e C°°(Q) with x = in \r\ < ^ and ^ = 1 in |r| > ^. If 
we let 

Vn,k - = X u n,k' 
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for < k < cin, and extend by outside of the hat of Q, then v n ^ are quasimodes for the Dirichlet 
Laplacian on Q with 

(-A-a^K^ = 0(2-" r/3 ). 

— C v /3 

Now, by the orthogonality of u n ^, v n ^ are orthogonal up to 0(e ). Hence, {v are a family 
of c 2"( c 2 n+1 ) a l mos t orthogonal quasimodes with 0(e~ Cn ' 13 ) error. Thus, 

{v n ,k' a l,k'° ^ k< c ^ n ' n = 

form a positive density set of quasimodes for -A on [0, 1]. 

Since the v n y are 0(n~°°) quasimodes, WF/,(o n ^) is invariant under the Hamiltonian flow (see, 
for example, fT7[ Section 12.3]). Therefore, since WF/,(t> n/ fc) does not intersect the foot of Q, the 
quasimodes must concentrate away from the ergodic set, U. Thus, v n ^\u has a defect measure y.. 
Hence, Theorem [3] applies and the constant we obtain is 0. 
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